Diffusion on non exactly decimable tree-like fractals 
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■ We calculate the spectral dijiHMi-iion of a. wide class of tree-like fractals by ,«>lviu,i; flic random 

walk problem through a new analytical technique, based on invariance under generalized cutting- 
decimation transformations. These fractals are generalizations of the NTd lattices and they are 
characterized by non integer spectral dimension equal or greater then 2, non anomalous diffusion 
f-H | laws, dynamical dimension splitting and absence of phase transitions for spin models. 
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I. INTRODUCTION 



The spectral dimension d of non translation invariant structures is up to now the best generalization of Euclidean 
dimension of lattices when dealing with dynamical and thermodynamical properties. It can be defined according to 
the large time asymptotic behavior of random walks and can be shown to be relevant for several different physical 
phenomena, such as vibrational dynamics, electrical conductivity and phase transitions. Many properties of d has 
been argued and proved starting from exact analytical calculations on a few particular cases and any further progress 
in understanding its relevance strongly depends on the availability of exact results in a wider range of particular cases. 
The most used mathematical tools to calculate d belong to two main classes: renormalization group and combinatorial 
techniques. Unfortunately renormalization group can give exact results only on exactly decimable fractals which, in 
turn, have been shown to have d < 2. On the other hand combinatorial techniques, based on iteration of cuttings, 
can be applied only to discrete structures with a given characteristic scale. Recently, a particular combination of both 
techniques has been successfully applied to the more complex case of NTd lattices Q . These lattices are not exactly 
decimable but they are invariant under a more complex geometrical transformation we shall call Cutting-Decimation. 
On NTd lattices the random walks problem has been analytically solved using a Cutting-Decimation method based on 
a time-rescaling technique H . The NTd has been shown to have remarkable properties such as non integer spectral 
dimension equal or greater then 2, non anomalous diffusion laws, absence of phase transitions || and dynamical 
dimension splitting Q . Because of these peculiar characteristics NTd lattices have been widely used for the study of 
statistical mechanics on non translation invariant lattices and they opened the way to the research of other structures 
sharing the same properties. 

In this paper we deal with generalizations of NTd lattices we shall call 2 m NTD, uNTd and p — polygonNTo 
introducing three corresponding new techniques to solve analytically the random walks problem. As a result we 
obtain very general structures with the same properties of NTd and, at the same time, we test new techniques that 
can be useful in other cases of non exactly decimable fractals. 



II. RANDOM WALKS ON 2 M NT D VIA EXACT TIME RESCALING 

The 2 m NTD are infinite fractal trees that can be recursively built using the following recipe (Fig. 1). An origin 
point O is connected to point A by a link of length 1 (the log of the tree); from A the log splits in k branches of 
length 2 m (i.e. made of 2 m consecutive links) which, in turn, split in k branches of length 2 2 ™ 1 and so on in such a 
way that each branch of length 2 nm splits in k branches of length 2( n+1 ) m . The simplest case is that of m — 1 and 
corresponds to the NTd lattices. The 2 m NTD invariance under the Cutting-Decimation transform can be described 
as follows (Fig. 2). Suppose to cut the log of the tree and to separate the k branches starting from A. Then each 
of these branches is exactly the same as the original tree with a dilatation factor equal to 2 m and therefore it can 
be reduced to the original tree through decimation. Let us now consider the discrete time random walks problem on 
2" 1 NTd, in order to calculate the spectral dimension d of the lattice. The spectral dimension of a graph is defined 
through the relation ||: 

Po(t)^t-^ 2 (1) 

where Po[t) is the probability for a random walker to return to the starting point O after a walk of t steps for t — > oo. 
The cutting transformation applied to random walks gives a relation between PQ Ce (t), the probability of returning to 
the starting site O after a t steps walk on the whole tree, and P^ ranch (i), the probability of returning to the starting 
site A after a t steps walks on one of the branches starting from A. This relation has been obtained in || for the case 
to = 1 in terms of the generating functions (discrete Laplace Transforms) P(X) of the probabilities P(t), but it holds 
for 2 m NTD for every to: 

- — branch 

Po tree (X) = PA b W h +k (2) 
(1 — A 2 )Pa (A) + k 

In the to = 1 case the decimation transformation has been performed using a time-rescaling technique. Indeed, the 
motion of the random walker on the branch considered only after an even number of steps can be exactly mapped in 
the motion of a random walker on the tree after the introduction of a staying probability pa — 1/2 in every site i. 
This equivalence can be translated in terms of generating functions through the substitutions : 

MX) - ^tPo (jtV) (3) 
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A ^ A 2 (4) 
In the case m = 1 equations (J3J) and ([|) can be used to rewrite (||) as: 



2 —tree f .2 \ 



tree , , 2-A 2 J ^ \ 2-A a / "T ^ 

(A) = - \ (5) 



If m > 1 this procedure must be iterated m times obtaining 

nr=i^)^^ ,i (A !+1 )+fc 



with 



— tree llltl 2^ , 

P ° (A) = V 7 ~ '\ X ~t.ee (6) 

(1 - A 2 ) (n™i i^A?) P o (Am) + fc 




i being the iteration step. To find the value of d we consider the singularities of the generating functions as A = 1 — e, 
e — > + . In terms of e the decimation transformation can be resumed as: 



as e — > 0, so that (||) becomes: 



- — branch - — tree 

P A (e)~2 m P (2 2m e) (7) 



2e2 m P (2 2m e) + jfc 

From (^), applying standard Tauberian theorems techniques ^] one obtains for a 2 m NT£> graph: 

~ In fc 

d 2m = 1 + — — (9) 
ln2 m w 

which represent the generalization of the result obtained for m — 1 . 

III. RANDOM WALKS ON NNT D AND P - POLYGON NT D VIA ASYMPTOTIC DECIMATION 

The previous results can be extended to nNTjj, where now n is an integer and not necessarily a power of 2, and 
to p—polygonNTu, where the branches of NTr> are replaced by p- vertices regular polygons (Fig. 3). Let us consider 
nNTu first. While relation ^ still holds, the exact time-rescaling procedure can not be applied to the branch of 
generic length n. However even in this case it is possible to obtain an asymptotic recursion relation applying the 
Renormalization Group techniques usually implemented on exactly decimable fractals (6). Although this procedure 

— - tree 

cannot give an exact equation for Po (A) as in the previous case, nevertheless it can be used to obtain the exact 
value of d via an asymptotic estimation. Indeed, in this case the branch of the nNTo can be considered as a tree with 
a dilatation factor equal to n. The log of this tree can be reduced to a unitary length log after the suppression of the 
n — 2 sites between the edges and introducing a new link connecting the edges. The same operation can be repeated 
for branches of every length suppressing the inner n — 2 consecutive sites in every sequence of n sites and introducing 
a new link between the surviving points. The final structure is equal to the original tree and the generating function 

branch — - branch 

Pa (e) becomes Pa (e ) where: 

e' = n 2 e (10) 

'branch , 1 — branch 

P A (e') = -PA (e) (11) 
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- — 'branch _ _ - — tree _ i— i 

Now Pa (e ) coincides with Po (e') since our branch has been transformed into a tree and (|2j) can be rewritten 

as: 

Po tree (e) = nF Cef £) + k (12) 
2enP (n 2 e) + k 

Using the procedure described in the previous section for 2 m NT[), from ( |l2| ) it follows that for an n — NTjj the 
spectral dimension is given by: 

~ In fc , . 

d n = l + - 13 

Inn 

An analogous technique can be used for p—polygonNTo (Fig. 3). The log polygon has now p faces of unitary length; 
from each of p — 1 of its vertices k polygons depart, whose faces have length n and so on. These structures, though 
similar to NTjj are no longer loopless structures nor necessarily bipartite graphs (e.g. the 3— polygon tree). The 
Cutting-Decimation transform can be applied to p—polygonNTp as in the case of NTjj with the same substitutions 
@ and (@. Indeed, even if (||) does not hold in this case, a new relation between the generating functions of the 
tree and that of one of its branches can be obtained using bundled structures theory [|7|. A bundle structure is 
composed graph obtained by joining to each point of a "base graph" a copy of a "fibre" graph in such a way that every 
fibre has only one point in common with the base and no points in common with the other fibers. Let us consider a 
p—polygonNTp and suppose to attach k branches also in the free vertex of the log (the root of the tree): we obtain 
a bundled structure having the log polygon as base and the graph made of k branches as fibre. Since for a p-polygon: 



a 



PoW ~ ^ d4) 



as A — + 1, we obtain for our bundled structure: 



* rw %_ A g-~'J l 1 -tTI .---». ' (15) 

— b.s. 

where Po (A) is the generating function of the probability of returning to point O (one of the vertices of the log 

— -branch 

polygon) after a random walk on the bundled structure and Fa (A) is the generating function of the probability 
of returning for the first time to the point of connection with the base after a random walk on the fibre. For a generic 
structure the generating function of the probability of returning to the starting point O, -Po(A) and of returning for 
the first time to the starting point, Fo{X), satisfy B: 

Po{\)= [l-Foi^y 1 (16) 

Using: 

b.s., , k ' — branch, , 1 ~ tree , , 

F ° {x) = kTT FA w + kTi Fo (A) (17) 

t A — I I — K A — — " t ree — ' branch 

where Fo (A) refers to the p-polygonNT D , from (|15|), ( |16[ ) and ( J17| j a relation between P (A) and Pa (A) 
follows, which represents the cutting transformation. It is now possible to perform the Cutting-Decimation transform 
to p—polygonNTn and get: 

~ In k(p — I) , 

d p = 1 + ^ '- 18 

Inn 

In the same way we can calculate the spectral dimension of an NTd built with d— dimensional simplexes instead of 
p— polygons. A G?-dimensional simplex is a complete graph of d + 1 points i.e. a graph where each point is nearest 
neighbor of all other points. The 2— dimensional case is the triangle, the 3— dimensional one is the tetrahedron and 
so on. Since for d— simplex Po(ty ~ l/(d + 1)(1 — A) the spectral dimension is: 

~ In kd , , 

d d = 1 + (19) 

Inn 



1-^ W P \ fc + i l-^ (A), 
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IV. CONCLUSIONS 



The generalized NTjj lattices described here have in general non integer spectral dimension depending on the 
geometrical features of the lattices, such as the growing factor 2 m or n, the number of branches k and the number 
of vertices of the polygon. Since the intrinsic fractal dimension always coincides with the spectral dimension, the 
diffusion ]2| on all these lattices is not anomalous i.e. is described by the asymptotic law: 

<r 2 {t)>~t a (20) 

with a = 1. Moreover as in the case of simple NTp lattices [Q, we can show that also generalized NTd present 
dynamical dimension splitting. This means that the vibrational spectral dimension d, which characterizes the density 
of vibrational modes as uj — > through the relation: 

p(w) ~ uj 1 - 1 (21) 

is different from the diffusive spectral dimension of random walks d. In particular we have d = 1 and this can be 
intuitively understood noting that the topology of generalized NTo is dominated by linear chains which become longer 
and longer in the outer branches. Since d = 1, by the generalized Mermin and Wagner Theorem |^], phase transitions 
with spontaneous breaking of a continuous symmetry for non zero temperature cannot occur on these structures. All 
these remarkable properties, typical of simple NTrj lattices, have been here obtained even in absence of some peculiar 
characteristics of NTry such as the bipartite and loopless nature of the graph. This enlarged family of lattices can be 
used, as NTjj lattices, to reach a better understanding of the geometrical features giving rise to dynamical dimension 
splitting and, at the same time, they represent the testing ground to study the relation between spectral dimensions 
and physical phenomena. 
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Figure captions 



Fig.l 2 rn NT D with k = 3, m = 1 

Fig. 2 Cutting-Decimation procedure: 

a) Cutting of the log of the NT D 

b) Separation of the k branches 

c) Decimation of the points labelled by X 

d) Recovering of the original NTrj 

Fig. 3 4 - polygonNT D with k = 1, n = 2 
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